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Abstract 

Using the matrix product formalism, we define a multi-parameter family of 
spin models on one dimensional chains, with nearest and next-nearest neighbor 
anti-ferromagnetic interaction for which exact analytical expressions can be found 
for its doubly degenerate ground states. The family of Hamiltonians which we 
define, depend on 5 continuous parameters and the Majumdar- Ghosh model is 
a particular point in this parameter space. Like the Majumdar- Ghosh model, 
the doubly degenerate ground states of our models have a very simple structure, 
they are the product of entangled (dimer) states on adjacent sites. In each of 
these states there is a non-zero staggered magnetization, which vanishes when 
we take their translation-invariant combination as the new ground states. At the 
Majumdar- Ghosh point, these entangled states become the spin-singlets pertain- 
ing to this model. We will also calculate in closed form the two point correlation 
functions, both for finite size of the chain and in the thermodynamic limit. 
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1 Introduction 



One of the basic problems of many body physics is to use analytical and numerical 
tools to find the ground state of a many-body system governed by a given Hamilto- 
nian. Unfortunately, except for a few exactly solvable examples, the task of finding 
the exact ground state of a given Hamiltonian is notoriously difficult. Moreover in 
some cases also like the anti- ferromagnetic Heisenberg chain, although the ground 
state is known, its structure turns out to be quite complicated, making the calcula- 
tion of correlation functions very difficult. 

In dealing with difficult problems, one way round the difficulty is to investigate 
the inverse problem, that is to start from states with pre-determined properties and 
find the family of Hamiltonians for which such a state is an exact ground state. The 
suitable formalism for following this path is the matrix product formalism, which was 
first realized in the work of [1] , then formalized in finitely correlated states [21 [3] , also 
known as optimal [U [5] or matrix product states. This formalism has been followed 
in recent years both constructing various models of interacting spins on chains, [HE], 
ladders [6l [7] and some two dimensional lattices [8] . 

It has also been studied to gain a better understanding of many phenomena in corre- 
lated systems, i.e. quantum phase transitions [9], renormalization group algorithms 
[lOj . and entanglement [llj . 

Of particular interest to us in this paper is the so called Majumdar- Ghosh model 
defined on a ring of 2N sites by the following Hamiltonian 

2N 



which has been introduced and solved exactly in a series of papers long ago [121 [131 HI], 
before the advent of the matrix product formalism, and later found to have a matrix 
product representation [15]. It was shown in the original papers that the ground 
state of this system is doubly degenerate, and each of these states has a very simple 
structure, that is, it is a product of singlets on adjacent sites. Let us call these 
ground states |(/>i) and \4>2)-, and denote a singlet state on two adjacent sites by 
S = "^(|01) ~ |10)), where |0) and |1) denote the up and down spin in the z direction, 
i.e. |0) = 1^:,+), |1) = 1^,— ). Then the degenerate ground states of ([T|) are: 



with -f^o|0i,2) = —^N\(t>i,2) ■ From these two states, one can form the new states 
l^±) '■= =t \4>2)^ which are eigenstates of the single- unit translation operator T 
on the ring, i.e. T\^±) = ib|^'-|-). The correlation functions on these new translation 
invariant ground states, denoted by the subscripts it, are defined as: 




(1) 



i=l 



\4>i) 

\^2) 



5'l,2 5'3,4 ••• -S'2Af-l,2Ar, 
5'2Af,l'S'2,3 • • • S2N-2,2N-l 



(2) 




a = X, y, z. 
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These were calculated in thermodynamic limit (N — > oo) in [12] and were found to 
be ^ 

(o-iOi = -2^r,2 a = x, y, z. (3) 

Since then many works have been done on different aspects of this modelfTG], including 
studies of spontaneous dimereization, excitation spectrum [17] and renormalization 
group study [IB] of this or very closely related models (i.e. models with anisotropy). 
It has also been shown [19] that the Majumdar- Ghosh model is the first member of a 
discrete series of models, all having exact dimer ground states, where exchange cou- 
pling decreases linearly with the separation of spins, and vanishes beyond a certain 
range, v (ly = 2 for the Majumdar- Ghosh model). 

The aim of this paper is to study the Majumdar Ghosh model in the context of 
matrix product formalism and in this way introduce a a multi-parameter deformations 
of this model and prove that its ground state is still doubly degenerate and dimerized. 
As is clear from ([T|), the Majumdar-Ghosh model is rotationally invariant. What we 
will do in this paper is that we reduce this symmetry to rotations in the x — y plane 
only and ask what are the Hamiltonians on the chain with nearest and next-nearest 
neighbor interactions, which are still exactly solvable. 

We will find a multi-parameter family of such exactly solvable Hamiltonians and 
show that their ground states are doubly degenerate and of dimer type. Both the 
Hamiltonian and their ground states will reduce to those of the Majumdar-Ghosh 
model in a certain limit. We will calculate all the correlation functions not only in 
thermodynamic limit, but for any finite A^. Our generalization of this model may shed 
light on some aspects of the Majumdar-Ghosh model and provide different proofs for 
some of its properties. 

Moreover, since the model is less symmetric and has 5 tunable parameters, it may be 
more flexible for matching with real experimental situations. 

The structure of this paper is as follows: first we review briefly the matrix product 
formalism in section ^ and then in section ([3]) we derive the Majumdar-Ghosh model 
by using the matrix product formalism. In section ^ we use the same formalism to 
construct a continuous multi-parameter family of models, and show that the original 
Majumdar-Ghosh model corresponds to a specific point in this parameter space. In 
various subsections of this section we make a detailed study of this model, derive its 
Hamiltonian, determine the structure of its ground state and prove that the doubly 
degenerate ground states of this model have a simple structure, that is, they are prod- 
ucts of entangled spin states on adjacent sites, in the same way that the ground states 
of the Majumdar-Ghosh model are products of singlets on adjacent sites. We also 
calculate all the one and the two point- functions of this model. In particular we show 
that each of the degenerate ground states has a non-zero staggered magnetization, 
however the translation-invariant ground states have zero staggered magnetization. 
We conclude the paper with a discussion. 
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2 Matrix Product States 



Let us first make a quick review of tlie matrix product states in a language wliich we 
find convenient [U]. For more detailed reviews of the subject, the reader can consult 
more comprehensive review articles [15]. Consider a ring of sites, where each site 
describes a d— level state. The Hilbert space of each site is spanned by the basis 
vectors i = 0, - ■ ■ ,d — 1. A state 

is called a matrix product state if there exist matrices Ai, i = 0, • • • , d — 1 (of 
dimension D), such that 

where Z is a normalization constant. This constant is given by 

Z = tr{E^), 

where 

d-l 
1=0 

Note that we are here considering homogeneous matrix product states (MPS) where 
the matrices depend on the value of the spin at each site and not on the site itself. 
More general MPS's can be defined where the matrices depend also on the sites [15j . 
It has been shown that such states can represent any quantum state on a chain [?] . 

The collection of matrices {Ai} and {fiUAiU~^}, where fi is an arbitrary complex 
number, both lead to the same matrix product state, the freedom in scaling with /i, is 
due to its cancelation with Z in the denominator of This freedom will be useful 
when we discuss symmetries. 

Symmetries: The state |^') is reflection symmetric if there exists a matrix 11 and 
a scalar a such that 

Af = aUAiU-^ V i. 

It is invariant under spin flip transformation if there exists a matrix X and a scalar 
A, such that 

A^ = XXAiX-^ V i. 

Finally, it is time-reversal invariant if there exists a matrix V, and a scalar r such 
that 

A* = rVAiV'^, V i. 

In this article, we work with real matrices so that the last condition is automatically 
satisfied. Let us now turn to continuous symmetries. 
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Consider a local symmetry operator R acting on a site as R\i) = Rji\j) where 
summation convention is being used, i? is a d dimensional unitary representation of 
the symmetry. A global symmetry operator TZ := R^^ will then change this state to 
another matrix product state 

where 

:= RijAj. 

The state 1^*) is invariant under this symmetry if there exist an operator U{R) and a 
scalar A_r such that 

R^JAJ = XRU-HR)AiU{R). (5) 
Repeating this transformation puts the constraint 

Thus U{R) \s a. D dimensional projective representation of the symmetry R. In case 
that i? is a continuous symmetry with generators T^, equation ([5|), leads to 

{Ta\jAj = [A,,Ta]+ tXaAi, (6) 

where Ta and 7^ are the d— and £)— dimensional representations of the Lie algebra 
of the symmetry, and jJLa are numbers. In this paper we will restrict our study to the 
case of = 0. Equations ([5]) and ([6]) will be our guiding lines in defining states with 
prescribed symmetries. 

The Hamiltonian: Given a matrix product state, the reduced density matrix of 
one site is given by 

_ tr{{A*®Aj)E^-^) 
~ tr{E^) 

This density matrix has at least d — zero eigenvalues. To see this, suppose that 
we can find complex numbers Cj such that 

Ec.-4.=0- (7) 

j=0 

This is a system of equations for d unknowns which has at least d—D^ independent 
solutions. Any solution Cj gives a null eigenvector of pij. Thus the null space of the 
reduced density matrix of one site comprises a subspace of at least d — dimension, 
spanned by the independent solutions of equation ([7]). 

The same reasoning applies to reduced density matrices of k consecutive sites which 
are given by 

_ tT{[Al...Al®A,,...A,^)E^-^) 
Pii-ikji-jk tr{E^) 
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In this case the null space of the reduced density matrix of k sites is spanned by the 
solutions of 

d-1 

ji,---,jk=0 

The number of independent solutions of this system of equation is now df'-D'^. Thus 
for the density matrix of k sites to have a null space it is sufficient (but not necessary 
and this is a crucial point) that the following inequality holds 

Let the null space of the reduced density matrix be spanned by the orthogonal vectors 
{a = 1, ■ ■ ■ , s > — D^). Then we can construct the local hamiltonian acting 
on k consecutive sites as 

s 

h Xa\ea){ea\, 

a=l 

where Aq,'s are positive constants. These constants together with the parameters of 
the vectors \ea) inherited from those of the original matrices Ai, determine the total 
number of coupling constants of the Hamiltonian. If we call the embedding of this 
local Hamiltonian into the sites / to / + A; by then the full Hamiltonian on the 

chain is written as 

N 
1=1 

The state l^*) is then a ground state of this hamiltonian with vanishing energy. The 
reason is as follows: 

N 

mH\^) = tr{H\^)m =J2tr{hi,i+kPi,i+k) = 0, 

1=1 

where pi^k+i is the reduced density matrix of sites / to / + A; and in the last line we have 
used the fact that h is constructed from the null eigenvectors of p for k consecutive 
sites. Given that H is a positive operator, this proves the assertion. 

Correlation functions: Let O be any local operator acting on a single site. Then 
we can obtain the one point function on site k of the chain (^'|0(A;)|^) as follows: 

(*|OWI*) = "•'^t.y"' . (9) 

where 

d-1 

Eo:= J2ii\0\j)A*^A,. 

i,j=0 

In the thermodynamic limit N — > oo, equation ([9]) gives 

{>^max\Eo\ 



mo\^) 



5 



where we have used the translation invariance of the model and Xmax is the eigenvalue 
of E with the largest absolute value, for which the matrix element is non-vanishing, 
and \\max) and {Xmax\ are the right and left eigenvectors corresponding to this eigen- 
value, normalized such that {Xmax\Xmax) = 1- Here we are assuming that the largest 
eigenvalue of E is non-degenerate. 

The n-point functions can be obtained in a similar way. For example, the two-point 
function {^\0{k)0{l)\^) can be obtained as 

(*io(.)o«i*) = "-'^"i^^ifg;"^"' (10) 

where Eo{k) := E^-^EqE'^. Note that this is a formal notation which allows us 
to write the n-point functions in a uniform way, it does not require that E is an 
invertible matrix. In the thermodynamic limit the two point function turns out to be 

mO{l)0{r)\^) = -^Y.^T^^>^^--\Eo\Xi){K\Eo\\max). 
■^max j 



3 The Majumdar- Ghosh model 

We take the size of the ring an even number 2N. Let us also choose d = 2 (spin 1/2 
particles), k = 2> (interactions up to next-nearest neighbors) and work with D = 3 
dimensional matrices. Note that the condition df^ > is not satisfied here, but in 
view of the fact that this is a sufficient and not necessary condition, this will not be a 
source of concern. The basic symmetry that we demand is the rotational symmetry 
around the z axis, which according to ([6]) demands that there be a matrix Sz such 
that 

[Ao,Sz] = ^Ao, [A,,Sz] = -^A,. 

Working in a basis where Sz = diagonal{— 1/2,0, 1/2), these equations immediately 
show that Aq and Ai should be of the following form: 




A, 




By a similarity transformation Ai — > SAiS~^ where S = diagonal{l,p,pq), the 
parameters of the matrix Aq are both set equal to 1, (the parameters of Ai change 
which we again denote them by g and h.) At this stage no similarity transformation or 
re-scaling, can reduce the number of parameters, however we note that the matrices 
Aq and Ai act as raising and lowering operators on the 3 dimensional space they act 
on, (z.e.^oN) = 1^ + 1)) Ai\i) oc |i — 1)) and hence the trace of any string of operators 
Ai-^Ai^ ■ ■ ■ Ai^j^ is non-zero only if the number of and Ai operators are equal in 
such a string. This implies first that any such string should be of even length (that 
is why we have taken a ring of even length), and second that a factor of /i^ will be 
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removed from the numerator and denominator of due to the normalization of the 
state. 

Thus the final matrices which will make a matrix product state with Sz symmetry, 
will be of the form: 




A, 
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This is all we can go with Sz symmetry and we will in fact make a detailed study of 
such models in the next section. 



What happens if we impose the discrete spin-flip symmetry on the state? This 
demands that there be a matrix X and a scalar A such that 



Aq = \XAiX~^, 



Ai = XXAqX-^. 



A short calculation shows that the matrix X which satisfies the above equation is 
equal to 

/ 1 



X 



A 

v2 



\ A^ 

where = g = 1 . Thus spin-flip symmetry requires that the parameter be a 
discrete parameter which we rename as g = e, where e = ±1. 



It is now found that the parity symmetry imposes no new condition on the pa- 
rameters, since with the following matrix 



n 



one can verify that = y/eliAiUr^ . Therefore the final form of the matrices leading 
to a model with a continuous Sz symmetry plus two discrete parity and spin flip 
symmetries will be of the form: 




Ao:-- 




A, 




[12) 



How is the matrix product state (jj]) with the above matrices is related to the 
states 101.2)? We will show in the next section, when we study the generalization 
of this model, that the matrix product state is equal to + |</'2)) that is it is the 
linear superposition of the ground states which is invariant under translation (with 
eigenvalue 1). In fact it has been shown that both \(j)i) and |(/>2) are ground states 
and so translation invariant combinations of them, namely \^±) := \4>i) ± \4>2), are 
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also ground states. The point is that only |^'+) has a matrix product representation. 
It sometimes happens that from a single matrix product state, one can derive the 
degenerate grounds states of a model not all of which have matrix product represen- 
tations. The Majumdar-Ghosh model is one such case. 

Let us also find the Hamiltonian. With the matrices (jl2p we should solve the 
equations X^i j,fc=o CijkAiAjAf^ = 0. It is readily found that the solution space of this 
system is spanned by the following four vectors: 

|ei) = |000) 

|e2) = |111) 

leg) = |001) - e|010) + |100) 

|e4) = |011) - ejlOl) + jllO). 

These vectors are parity invariant (invariant under left-right reflection) and are mapped 
to each other under spin- flip transformation (|0) — > |1), |1) ^ |0)). A Hamiltonian 
with all three symmetries is constructed as: 

h = J(|ei)(ei| + |e2)(e2|) +-fC(|e3)(e3| + |e4)(e4|), 

which when expanded in terms of the spin operators, and embedded into the chain, 
makes the total Hamiltonian equal to the following: 

^ K J+{2e-l)K J-m 
-n = 2^ --f^eo-j • CTj+i + —(Ji ■ crj+2 H ^ o-i,zO'i+i,z H ^ o'i,20'i+2,^. 

i=l 

For e = 1 this reduces to 

rr ^ ^ J + ^ J - 3K 

-ni = 2^ --fVfJj • CTj+i + —ai ■ (Tj+2 H — CTj.zCTj+i^^ H ai^zai+2,z, 

i=l 

and for e = — 1 reduces to 

= 2^J^<^i- <^i+l + -^Cn ■ H ^ {2cri,zO'i+l,z + <yi,z<yi+2,z)- 

i=l 

It is the Hamiltonian H^i which turns into the already familiar Majumdar-Ghosh 
model by taking J = S-fT. For general values of J and K it describes the 
anisotropic Majumdar-Ghosh model, where the anisotropy parameter j '^j is the 
same for both the nearest and next-nearest neighbors. 

We do not know whether this is already well-known that this model is also exactly 
solvable. In fact one can check directly that the eigenstates ([2]) are also eigenstates 
of the operator 

2N 

Hz ■■= 51 ^CTj^CTi+l + '7i(^i+2^ 
1=1 
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with eigenvalue 2N . To see this one needs to invoke the following relations which 
hold true for any two consecutive singlets 5i2 and 534: 



^1 '^2*^1,2 5*3,4 — '^la^iA 



'^2 '^3*^1,2 •5*3,4 = "^^1,2^^3,4) 



where U = ^(|01) + |10)). 

Using these relations one sees that in -^2 101,2) all the terms containing the U terms 
cancel out and one is left with Hz\(j)i^2) = —'2^\4'i,2)- 

Putting all these together we find 

— T — "^K 

H-i\ct>i,2) = ^ A^|<^1,2). 

What is the Hamiltonian Hil Does it represent any new model? The answer is 
negative. To see this we note that by rotating the even-numbered spins around the z 
axis by an angle of vr, under which 02]^ — > ~'^2i^ and cifj — > crfj, while keeping the 
odd-numbered spins intact, the Hamiltonian Hi is mapped onto H_i. 
This transformation will map the states |(/>i,2) onto the states \xi,2) where 



IXl) — f^l, 2^^3,4 ••• U2N-1,2N, 
1X2) = U2N,lU2,3- ■ ■U2N-2,2N-l- 

These states are the doubly degenerate ground states of Hi. 

4 Deformation of the Majumdar-Ghosh model 

In the last section we have shown how the Majumdar-Ghosh Hamiltonian and its 
ground states are derived in the matrix product formalism. Let us now relax the con- 
dition of parity symmetry, which automatically removes also the spin flip symmetry. 
Then the matrices will be of the form (jllh . which still have the continuous symmetry 
of rotation around the z axis. It is now easily verified that the solution space of the 
system of equations ([8]) with /c = 3 is spanned by the following four vectors: 

1 000) 
1111) 

|001) -5IOIO) + 5^1100) 
|011)-5|101)+5'|110). 

The most general local Hamiltonian having the matrix product state as its ground 
state is now given by 

h = Ji|ei)(ei| + J2|e2)(e2| + J3|e3)(e3| + J4\e4,){ei\. 



62) = 
es) = 
64) = 
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Rewriting this in terms of spin operators we find after some rather lengthy calculations 
the following form of the final Hamiltonian 

1 ^ 

H = - ^ Ki ai ■ (Tj+i + K2 ai- ai+2 + '^i'^i+i + -^40"f + K^o't + K^afafj^i 

i=l 

where the couplings are 





= -2g{l+g^){.h + J^) 


K2 


= 2g\.h + Ji) 


K3 


= 2(Ji + J2)+2(J3 + J4)(5-/+5' 


K4 


= Jl + J2-(J3 + J4)(l+/+/) 


K5 


= 3(Ji- J2) + (J3- J4)(l + 5' + /) 




= Ji-J2-{J-i-JA){i-9 + g^f 


K7 


= -25J3 + 25V4 


Ks 


= -25V3 + 25J4 


Kg 





This is a 5-parameter family of models with nearest neighbor and next-nearest neigh- 
bor interactions, whose ground state can be written down exactly. We will show that 
with the couplings as above, this model has a doubly degenerate ground state, each 
of which has a simple product structure of dimers, like the one given in ([2]). 

Note that the Hamiltonian (jl3p will reduce to H-.i when we restrict the parameters 
as follows 

g = -1, Ji = J2 = J, J3 = J4 = K 
and to ([1]) when we further set J = 3K. 



4.1 Structure of the matrix product state 

The matrix product state (jlip has a simple structure. We prove the following theorem: 

Theorem: The matrix product state ^ constructed from the matrices (llip which 
is a ground state of has the following structure. Define a state 

S{g) :=5|01) + |10), 
IMa)) = S^^{g)S^^{g)S.^^{g)--- S 

2JV-2,2JV-1 (5), (14) 

that is \(p2i9)) = T\(j)i{g)), where T is the translation operator along the chain by one 
lattice unit. Then we have 

\^{g)) = \M9)) + \M9))- 
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Proof: Consider the matrix product state \^{g)) = J2 ^jii2---i2jv Ni^2 • • • ^2Af)- The 
amphtudes ^'iij2---«2jv non-vanishing only when 

i) the number of indices and 1 indices are equal (in view of the fact that the ma- 
trices Aq and Ai are raising and lowering operators), ii) no three consecutive indices 
are or 1, (since Aq = Af = 0), and iii) no substrings of the types • • • 11011 • • • or 
• • • 00100 • • • appear in the indices, (since AjAoAl = A^AiA^ = 0.) 

Let us define the compact notations := 01 and 1 := 10. In view of the above 
three properties, we can divide the set of states ^2, • " " ^2Af) into two types, type 
1, which are states of the form |00110) = |0101101001), that is, strings of O's and 
l,s in arbitrary order, and type 2, which are transformed to the above type by a 
one-unit cyclic translation like 1 1010110100). There is a one to one correspondence 
between these two types of states and hence the state \^) can be divided in two parts, 
depending on the type of basis states {"^{g)) = |(/)i(g)) + |(/>2(fi')) , where 

<^l(5) = X! </'aia2-aivl«l«2 • • • ttAf) 

ai,«2,---:"iV=0,l 

and \(j)2ig)) = T\(l)i{g)). We now note that 

Ao = AoAi =1 g \ Ar = A^A^ =\ 1 | . (15) 

This means that ^-nd A-i commute with each other, which immediately implies 
that 101(5')) is a product state, i.e. 

N N 
\<Pi{g)) = E tr-(^-Af-'")|0-l^— ) = E 5™|0™1^— ) = (5|0) + |1))^^, 

m=0 m=0 

where lo^l^"*") denotes a uniform superposition of states each of which has m O's 
and N — m I's, in different positions. This completes the proof. 



Are the two individual states \(t)i{g)) and \4)2{g)) the ground states of HI There 
are supports for thinking so. First we have checked this for low values of system 
size and second we know that as g — > —1, these states approach the states 
of Majumdar-Ghosh ^ which are known to be ground states of the Hamiltonian 
H-i. Therefore we can make two linear combinations of these ground states, namely 
\^+{g)) := |(?!)i(5)) + |(/)2((7)) which has a matrix product representation and |^_((7)) := 
\'t'i{g)) ~ \4'2{g)) which does not have such a representation. In the next subsection 
we calculate the correlation functions on these translation invariant states and they 
again approach the ones calculated in [T2] in the limit g = —1 and N — > oo. For our 
future reference we need the following easily verified properties: 

{M9)\M9)) = {M9)\M9)) = (1 + 5')"^ {<Pi{g)\M9)) = {M9)\M9)) = 25^. 
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4.2 Correlation functions 



In the Majumdar- Ghosh model both the ground states ^ have anti-ferromagnetic 
order and {cr^'^'^) = 0. This is clear from the fact that the states ([2]) are product of 
singlets on adjacent sites. In our model, the states are no longer singlets and while 
there is no magnetic order present, the staggered magnetization is non- vanishing. 
Consider the states \(j)i{g)), and \(j)2ig))- In view of their simple product structure, it 
is readily seen that 

(</.i(<7)k!l</>i(<7)) = -{M9)W2\M9)) = {Si2Wt\S^2) = (16) 



and 



{Mg)Wt\M9)) = -{Mg)WI\Mg)) = {S2N,iWt\S2N,i) = (i7) 



g 

t 

i + g' 



This means that the value of staggered magnetization defined as 



1 ^, 



i=l 

is non-zero in these two states, i.e. 

However there is no staggered magnetization in the translation invariant states 
l^'-t). To see this we can calculate af for the state |^'+) = {(j^iig)) + \4'2{9)), either 
by using the matrix product formula Q or by a direct calculation in which we use 



([HI), (HZD and 



in which O is any operator and we have abbreviated \(pi{g)) to \(pi) and used the 
fact that = {(p2\4>2) to simplify terms in the denominator. For this direct 

calculation, the only new quantity which is needed is {(j)i{g)\ai\(j)2ig)) ■ We know 
that all the states in the expansion of \(j)i{g)) and \<j)2{g)) are different from each 
other except the states |010---1) and |101---0) which are common to both and 
transformed into each other by the action of the translation operator T. So we find 

{M9)K\Ug)) = (5"(oi • • • i| + (10 • • • o|)|af |(5"|io • • • 0) + |oi • • • i)) = o, (19) 

implying that there is no staggered magnetization in \^±). 

The two point functions can also be calculated along similar lines. We point out 
only the general method. The product nature of {cpi) and \(j)2), formula ([TH]) and 
the reasoning used above for the calculation of the cross term {(t>i{g)\cria^\(p2{g)) 
facilitate calculations of {o'fa^)±. However the cross term (0i(5')|o-f (T^|02(fl')) is hard 
to calculate in this way. Instead we calculate {afa^)-^ (i.e. for the matrix product 
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state 1^+)) from (jlOp and infer the cross term from its result. We then use this cross 
term for determining the correlation {afcrf)-. The reader can verify these results 
by direct calculation for chains of small size. The interesting point is that all the 
correlation functions can be expressed in terms of a variable u := 5 + | which is 
invariant under the interchange of g — > |. The results are: 

2 ±u^{2u-^-l) , z z X f ,,,-2±n^(4'u-2 - 1) 



(^f^f )± = ' ^^'^'>2^ = ^"^^ 2 ± 

and 



I „n n\ - ^ I „n „ 



2 ±n~' \-i"'->2/- ^X^pv ■ 

It can be verified that these correlation functions reduce to ^ at the point g = —1 
[u = —2), after we take the thermodynamic limit {N — > oo). 

5 Discussion 

We have constructed a 5-parameter family of spin-Hamiltonians with nearest and next 
nearest neighbor interactions and have found their exact ground states. The ground 
states, depend on only one of these parameters, are doubly degenerate and have a very 
simple structure. They are product of entangled states (dimers) on adjacent spins. We 
have also calculated the one and the two point functions in closed analytic form both 
for finite size of the chain and in the thermodynamic limit. The dimerized ground 
states show staggered magnetization, but the translation-invariant combinations of 
these states have no staggered magnetization. 
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